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ABSTRACT 

Using a simple off-axis jet model of GRBs, we can reproduce the observed 
unusual properties of the prompt emission of GRB 980425, such as the extremely 
low isotropic equivalent 7-ray energy, the low peak energy, the high fluence ratio, 
and the long spectral lag when the jet with the standard energy of ~ 10 51 ergs 
and the opening half-angle of 10° < A8 < 30° is seen from the off-axis viewing 
angle 9 V ~ A6 + IO7 -1 , where 7 is a Lorentz factor of the jet. For our adopted 
fiducial parameters, if the jet that caused GRB 980425 is viewed from the on-axis 
direction, the intrinsic peak energy E p (l + z) is ~2. 0-4.0 MeV, which corresponds 
to those of GRB 990123 and GRB 021004. We also discuss the connection of 
GRB 980425 in our model with the X-ray flash, and the origin of a class of GRBs 
with small Ej. 

Subject headings: gamma rays: bursts — gamma rays: theory 



1. INTRODUCTION 

Recently, a very luminous gamma-ray burst (GRB), GRB 030329 at the distance of 
0.8 Gpc (z = 0.1685) was confirmed to be associated with supernova SN 2003dh (Stanek 
et al. 2003; Uemura et al. 2003; Price et al. 2003; Hjorth et al. 2003). The geometrically 
corrected 7-ray energy E 1 of this event ~ 5 x 10 49 ergs is a factor 20 smaller than the standard 
value, if the jet break time of ~ 0.48 days is assumed (Vanderspek et al. 2003; Price et al. 
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2003). GRB 980425 was the first GRB associated with a supernova (SN) event, SN 1998bw 
at z = 0.0085 (36 Mpc) (Galama et al 1998; Kulkarni et al. 1998; Woosley et al. 1998; Pian 
et al 2000, 2003). There are some other events that might be associated with supernovae 
(Delia Valle et al. 2003; Wang & Wheeler 1998; Germany et al. 2000; Rigon et al. 2003). 
Therefore the association of the long duration GRBs with supernovae is strongly suggested 
and at least some GRBs arise from the collapse of a massive star. 

In this context, it is important to investigate whether GRB 980425/SN 1998bw is sim- 
ilar to more or less typical long duration GRBs like GRB 030329/SN 2003dh. However, 
GRB 980425 showed unusual observational properties. The isotropic equivalent 7-ray energy 
is Ei SO ~ 6 x 10 47 ergs and the geometrically corrected energy is E 1 ~ 3 x 10 46 ergs (A0/O.3) 2 , 
where A8 is the unknown jet opening half-angle. These energies are much smaller than the 
typical values of GRBs £ 7 ~ 1 x 10 51 ergs (Bloom et al. 2003; Frail et al. 2001). Bloom 
et al. (2003) claim that there should be some events with small E 1 such as GRB 980519 
and GRB 980326 and that GRB 980425 might be a member of this class. The other prop- 
erties of GRB 980425 are also unusual; the large low energy flux, the long spectral lag, the 
low variability, and the slowly decaying X-ray luminosity of its counterpart detected and 
monitored by BeppoSAX and by XMM-Newton (Frontera et al. 2000a; Norris, Marani, & 
Bonnell, 2000; Fenimore & Ramirez-Ruiz 2000; Pian et al. 2000, 2003). 

Previous works suggest that the above peculiar observed properties of GRB 980425 
might be explained if the standard jet is seen from the off-axis viewing angle (Ioka & Naka- 
mura 2001; Nakamura 1999; Nakamura 2001; see also Maeda et al. 2002; Iwamoto 1999; 
Dado, Dar, & De Riijula 2003; Dar & De Riijula 2000). Following this scenario, the rela- 
tivistic beaming effect reduces E iso and hence E 1 . The quantity E iso is roughly proportional 
to 5 2 ~ 3 for the typical observed spectrum, where 5 = [y(l—(3 cos(9 v — A#))] _1 is the Doppler 
factor and 9 V is the viewing angle (Yamazaki, Ioka, & Nakamura 2002). Since E iso is ~ 10 4 ~ 5 
times smaller than the standard value, 5 should be 20 ~ 10 2 times smaller than the usual 
value. Then the peak energy E v (oc 5) becomes 20 ~ 10 2 times smaller than on-axis E p , that 
is measured when the jet is seen from the on-axis viewing angle. However, the observed E p of 
GRB 980425 (~ 50 keV) is only a factor 4 or 5 smaller than the typical value of ~ 250 keV. 
Therefore, one might consider that GRB 980425 belongs to a different class of GRBs. 

It is well known that the distribution of E p is log-normal with the mean of (E p ) ~ 
250 keV (Preece et al. 2000). Ioka & Nakamura (2002) showed that if the distribution of 
intrinsic E p (i.e. E p (l + z)) is log-normal, the redshifted one is also log-normal under the 
assumption that the redshifts of the observed GRBs are random. Therefore, (E p (l + z)) ~ 
570 keV since the mean value of the measured redshifts is ~ 1.3 (Bloom et al. 2003). There 
are some GRBs with even higher intrinsic peak energy; for example, E p (l + z) ~ 2.0 MeV 
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for GRB 990123 (Amati et al. 2002) and E p (l + z) ~ 3.6 MeV for GRB 021004 (Barraud 
et al. 2003). Furthermore, Fig. 3 of Schaefer et al. (2003) shows that the highest value of 
E p (l + z) detected by BATSE is about 4 MeV. Since GRB 980425 is the nearest GRB, the 
redshift factor is not important. In this sense, the peak energy of GRB 980425 is at least 
a factor ~ 10 smaller than the usual one of ~ 570 keV. Suppose that the intrinsic E p of 
GRB 980425 is similar to that of GRB 990123 and GRB 021004 when the jet of GRB 980425 
is seen from the on-axis viewing angle. Then, the observed E p of GRB 980425 is ~ 10 2 times 
smaller than the intrinsic E p of GRB 990123 and GRB 021004. This is the reason why we 
incline to reconsider the off-axis jet model for GRB 980425. 

In this Letter, assuming the rather large on-axis E p , we reconsider the prompt emission 
of GRB 980425 using the simple model in Yamazaki, Ioka, & Nakamura (2002, 2003b) to 
reproduce its unusual observed quantities. In § 2, in order to extract the observational 
properties that should be compared with our theoretical model for prompt emission of the 
GRB, we analyze the BATSE data of GRB 980425. In § 3, we describe a simple jet model 
including the cosmological effect. We assume a uniform jet with a sharp edge. Numerical 
results are shown in § 4. Section 5 is devoted to discussions. Throughout this paper, we 
adopt the flat universe with Q M = 0.3, Q\ = 0.7 and h = 0.7. 

2. SPECTRAL ANALYSIS FOR PROMPT EMISSION OF GRB 980425 

USING BATSE DATA 

In our simple jet model of GRBs, the time-dependence of spectral indices is not treated, 
while it is known that the spectral parameters of GRB 980425 changed in time (Galama 
et al. 1998; Frontera et al. 2000a). Hence, we should discuss the time-averaged observed 
spectral properties of GRB 980425 before we apply our model to them. 

Using the BATSE data of GRB 980425, we analyze the spectrum within the time of 
Full- Width at Half-Maximum of the peak flux in the light curve of BATSE channel 2 (50-110 
keV). This time interval approximately corresponds to portions "B" and "C" in Frontera 
et al. (2000a), when most of photons arrived at the detector and the spectral shape was 
approximately constant with time. We fit the observed spectrum with smoothly broken 
power-law function given by Band et al. (1993), that is characterized by the energy at the 
spectral break E , and the low- and high-energy photon index a and (3, respectively. For 
the case of (3 < —2, the peak energy is derived as E p = (2 + a)E . The best-fit spectral 
parameters are 



a = 



-1.0 ±0.3 
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(3 = -2.1 ±0.1 

E p = 54.6 ± 20.9 keV. 

The reduced chi square is 1.10 for 31 degree of freedom. These results are consistent with 
those derived by the previous works (Frontera et al. 2000a; Galama et al. 1998). Although 
the photon indices are the typical values of GRBs, E p is lower than the typical values of 
GRBs (Preece et al. 2000). This spectral property is similar to one of the recently identified 
class of the X-ray flash (Kippen et al. 2002; Heise et al. 2001). 

The observed fluence of the entire emission between 20 to 2000 keV is S(20— 2000 keV) = 
(4.0 ± 0.74) x 10~ 6 erg cm -2 , so the isotropic equivalent 7-ray energy becomes E iso = (6.4 ± 
1.2) x 10 47 ergs. The fluence ratio is R s = S(20 - 50 keV)/,S(50 - 320 keV) = 0.34 ± 0.036. 
In the following sections, we reproduce the above results using our prompt emission model. 



3. MODEL OF PROMPT EMISSION OF GRBs 

We use a simple jet model of prompt emission of GRBs adopted in Yamazaki et al. 
(2003b), where the cosmological effect is included (see also Yamazaki et al. 2002, 2003a; 
Ioka & Nakamura 2001). We adopt an instantaneous emission of infinitesimally thin shell at 
t = t and r = r . Then the observed flux of a single pulse is given by 



A } dl [ 7 (l-/W(T))f ' [L) 

where, 1 — /3cos6(T) = (l + z)~ 1 (c/3/ro)(T — To) and Aq determines the normalization of the 
emissivity. Detailed derivation of Eq. (1) and the definition of A0(T) are found in Yamazaki 
et al. (2003b). In order to have a spectral shape similar to that derived by the previous 
section, we adopt the following form of the spectrum in the comoving frame, 

(z/>o) 1+ ° s exp(-i//i/ ) for ^'K < "b " Pb (2) 

(u'/u' ) 1+ ^(a B - (3 B Y B ~^ exp(/3 B - a B ) for i//i/ Q > a B - (3 B , U 



with a B = — 1 and f3 B = —2.1. Equations (1) and (2) are the basic equations to calculate the 
flux of a single pulse, which depends on the following parameters; 7, 71/Q, 9 V , A6, ro/cft^ 2 , 
z, and Aq. In the next section, the viewing angle 9 V and the jet opening half angle A# are 
mainly varied. The other parameters are fixed as follows; the quantity 7 is fixed as 7 = 100. 
The isotropic 7-ray energy is calculated as E iso = Atc(1 + z)^ 1 d 2 L S(20 — 2000 keV), where 
S(ui — v 2 ) is the observed fluence in the energy range hui-hu 2 keV. We fix the amplitude 
Aq so that the geometrically-corrected 7-ray energy E 1 = (A6) 2 E iso /2 be observationally 
preferred value when we see the jet from the on-axis viewing angle 9 V = 0. It is shown that 
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E 7 is tightly clustering about a standard energy £ 7 of ~ 10 51 ergs (Bloom et al. 2003; see 
also Frail et al. 2001; Panaitescu & Kumar 2002). Bloom et al. (2003) derived this energy 
as 

log£ 7 = log [1.15 x 10 51 (/i/0.7)~ 2 ergs] ± 0.07 , (3) 

so that £ 7 =(0.98-1.35) xlO 51 ergs, at the 1 a level while £ 7 =(0.51-2.57) x 10 51 ergs, at 
5 a level. Note that the smaller jet opening half-angle A9 corresponds to the larger A 
(Yamazaki et al. 2003b). 

Practical calculations show that when the jet with a B = — 1 and /3 B = —2.1 is seen from 
the on-axis viewing angle 9 V = 0, the observed peak energy becomes Ep v ~ ^ ~ 1.54 71^(1 + 
z)^ 1 , which is independent on A9 larger than ~ 7 _1 . In order to reproduce the observed 
quantities of GRB 980425, we adopt the value 7^ = 2600 keV, which yields E^ v=0 \l + 
z) ~ 4.0 MeV. For comparison, we consider another case of 7z/q = 1300 keV, which reads 
E { p v=Q \l + z) ~ 2.0 MeV. These values correspond to the intrinsic E p of GRB 021004 
and GRB 990123, respectively. Note here that in our jet model the quantities that will be 
calculated in the next section do not depend on r /c/37 2 ; for example, E iso oc A (r / 'cf3^ 2 ) 2 oc 
(r /c/37 2 ) since A oc (r /c/37 2 ) -2 . The value of r /c/37 2 will be determined when we discuss 
the spectral lag in § 5. 



4. ISOTROPIC ENERGY, PEAK ENERGY, AND FLUENCE RATIO 

We now calculate the isotropic equivalent 7-ray energy E iso as a function of 9 V and AO. 
Then, the peak energy E p and the fluence ratio R s = S(20 - 50 keV)/,S(50 - 320 keV) are 
computed for the set of AO and 9 V that reproduces the observed E iso of GRB 980425. 

For fixed A# and £ 7 , E iso is calculated as a function of the viewing angle 9 V . The result 
is shown in Fig. 1. When 9 V < A9, E iso is essentially constant, while for 9 V > A^, E iso 
is considerably smaller than the typical value of ~ io 51 ^ 53 ergs because of the relativistic 
beaming effect. In order to explain the observation, 9 V should be ~ 21° in the case of 
A# = 15°, while 9 V ~ 25° in the case of A# = 20°. This result does not depend on 7^ so 
much. 

The upper panels of Fig. 2 and 3 show for which E iso becomes equal to the observed 
values, as a function of A# in the case of 7z/q = 2600 keV and 1300 keV, respectively. Since 
the emissivity (oc A ) of the jet is small for large A9, the relativistic beaming effect should be 
weak for large A9. Therefore, the value of 9* — A9 is a decreasing function of A9. For such 9*, 

(0 =6*) (0 0* } 

we calculate the fluence ratio R* = R s v ~ and the peak energy E* = Ef v -° v >. The middle 
and the lower panels of Figs. 2 and 3 show the results. The quantity E* is proportional to 
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the Doppler factor 5 ~ [7(1 — (3 cos(#* — AO))]^ 1 . Therefore, when AO increases, 9* v — AO 
decreases so that E* increases. Since we fix spectral indices «b and /3b, i?* depends only on 
E*. Hence, if E* is large, the spectrum is hard and R* is small. For the fiducial parameters 
of 7z/q = 2600 keV, £ 7 = 1.15 x 10 51 ergs, and E iso = 6.4 x 10 47 ergs, AO should be between 
~ 18° and ~ 31°, and then #* ranges between ~ 24° and ~ 35° in order to reproduce the 
observed values of R s and E p . When £ 1 is varied from 0.51 x 10 51 to 2.6 x 10 51 ergs (at 5 a 
level), the allowed region with 20° < AO < 30° can exist even in the case of 7^ = 1300 keV. 

Note that 7 does not affect our results for observed R* and E* When 7 is large, 0* 
becomes small because the observed flux for fixed V becomes small due to stronger relativistic 
beaming effect. However, we can see that 7(0* — AO) remains almost unchanged even if 7 is 
varied. Then for fixed 71^, E* remains constant since E* oc v' 5 ~ 27^ [1 + (7(6** — AO)) 2 ]^ 1 . 
The quantity R* depends only on E* so that 7 does not affect the estimate of R*. 

5. DISCUSSION 

We considered the time-averaged emissions, which means that successive emissions from 
multiple subjets (or shells) are approximated by one spontaneous emission caused by a single 
jet (Yamazaki et al. 2002). We choose a B = -1, Pb = -2.1, 7 = 100, and 7^ = 2600 keV 
for the canonical set of parameters. As a result, when the jet of opening half- angle of 
AO ~ 10-30° is seen from the off axis viewing angle of 9 V ~ AO + 6°, observed quantities 
can be well explained. Derived 9 V and A# are consistent with those suggested in Nakamura 
2001, Nakamura 1999, and Maeda et al. 2002. We may also be able to explain observed low 
variability since only subjets at the edge of the cone contribute to the observed quantities 
(see the discussion of Yamazaki et al. 2002). If the jet is seen from an on- axis viewing angle 
(i.e. 6 V < A9), the intrinsic peak energy E P {1 + z) is ~ 4.0 MeV, which is almost the same 
as the highest one (Schaefer 2003; Amati et al. 2002; Barraud et al. 2003). 

As we have mentioned in § 3, E iso , E* and R* do not depend on the parameter r Q /j3c~f 2 . 
In order to estimate the value of ro//5c7 2 , we discuss the spectral lag of GRB 980425 (Ioka & 
Nakamura 2001). In our model, we can calculate the spectral lag AT, which is defined, for 
simplicity, as the difference of the peak time between BATSE energy channel 1 and 3. We 
obtain AT/(r /c^ 2 ) =0.97-1.34. Therefore, the observed value of AT = 3 s (Norris et al. 
2000) can be explained when r /cf3^ 2 = (2.2-3.1) sec, which is in the reasonable parameter 
range. 

The observed quantities of small E p and large fluence ratio R s (see also Frontera et al. 
2000a) are the typical values of the X-ray flash (Heise et al. 2001; Kippen et al. 2002; see 
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also Barraud et al. 2003; Arefiev, Priedhorsky & Borozdin 2003). The operational definition 
of the X-ray flash detected by BeppoSAX is a fast X-ray transient with duration less than 
~ 10 3 seconds which is detected by WFCs and not detected by the GRBM (Heise et al. 2001). 
If the distance to the source of GRB 980425 that has an opening half-angle of A9 = 20° 
were larger than ~ 86 Mpc, the observed flux in the 7-ray band would have been less than 
the limiting sensitivity of GRBM ~5x 10~ 7 erg cm~ 2 in 40-700 keV band (Band 2003), so 
that the event would have been detected as an X-ray flash. 

We might be able to explain the origin of a class with low i? 7 , pointed by Bloom et 
al. (2003). Let us consider the jet seen from a viewing angle 9 V ~ A9 + 7~\ where 7, 
is the Lorentz factor of a prompt 7-ray emitting shell. Due to the relativistic beaming 
effect, observed of such a jet becomes an order of magnitude smaller than the standard 
energy (see Fig. 1). At the same time, the observed peak energy E v is small because of the 
relativistic Doppler effect. In fact, the observed E v of GRB 980326 and GRB 981226 are 
~ 35 keV and ~ 60 keV, respectively (Amati et al. 2002; Frontera et al. 2000b). In our 
model the fraction of GRBs with low E 1 becomes 2/{ r jiA9) ~ 0.1 since the mean value of 
A9 ~ 0.2, while a few GRBs with low E 1 are observed in ~ 30 samples (Bloom et al. 2003). 
In later phase, the Lorentz factor of afterglow emitting shock 7/ is smaller than 7^, so that 
9 V < A9 + 7J 1 . Then, the observed properties of afterglow may be similar to the on-axis 
case 9 V <C A9; hence the observational estimation of the jet break time and the jet opening 
angle remains the same. 
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Fig. 1. — The isotropic equivalent 7-ray energy E iso is shown as a function of the viewing 
angle 9 V for a fixed jet opening half-angle A9. The source is located at z — 0.0085. The 
values of A6 are shown in parentheses. Solid lines correspond to the case of 7^ = 2600 keV, 
while dotted lines 7^ = 1300 keV. Other parameters are fixed as a B = — 1, Pb = —2.1, 
7 = 100, and £ 7 = 1.15 x 10 51 ergs. Horizontal dashed line represents the observed value of 
GRB 980425 E iso = 6.4 x 10 47 ergs. The value of E iso in the on-axis case 9 V < A9 is slightly 
smaller for 7Z/Q = 2600 keV than for 7Z/Q = 1300 keV. This is because the amplitude Aq is 
fixed so that we should observe constant E 1 from the source at z — 1, and the K-correction 
is larger for ^u' = 1300 keV than for 7Z/Q = 2600 keV. 
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Fig. 2. — The upper panel shows #* for which E iso is the observed value of GRB 980425, 

(9 =0*) 

while the middle and the lower panels represent the fluence ratio R* = R s v and the peak 

(Q =Q*) 

energy E* — E p v v , respectively. Solid lines correspond to the fiducial case of Ei SO = 6.4 x 
10 47 ergs and £ 



^ 7 — 1.15 x 10 51 ergs. The dotted lines represent regions where E iso becomes 
(6.4 ± 1.2) x 10 47 ergs when £ 7 is in 1 a and 5 a level around the fiducial value, respectively. 
Other parameters are fixed as a B = —1, Pb = —2.1, 7 = 100, and 7z/q = 2600 keV. 
The dot-dashed line in the upper panel represents 9* v = A9. Horizontal dashed lines in 
the middle and the lower panels represent the observational bounds R s = 0.42 ± 0.13 and 
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Fig. 3.— Same as Fig. 2 but for 7^ = 1300 keV. 



